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Abstract 

We show that the elastic charge of colloids in a nematic liquid crystal can be generated by 
the vector of external torque. The torque components play the role of two component charge 
(dyad) and give rise to the Coulomb-like potential, while their conservation law plays the role 
similar to that of Gauss' theorem in the electrostatics. The theory is applied to the colloid-surface 
interaction. A wall with homeotropic or planar director is shown to induce a repulsive 1/r^ force 
on the elastic dipole. The external torque, however, induces the elastic charge in this colloid and 
triggers switching to the repulsion. 
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I. INTRODUCTION 



Particles of a submicron and micron size immersed in a nematic liquid crystal (NLC) 
interact via the director field n which mediates the distortions induced by their surfaces 
P El El H El El [a El El Uni [n]. The new field of nematic colloidal systems, or nematic 
emulsions [lO], has gained a continuous growing interest over the past few years. The 
physics of these anisotropic colloidal systems has a deep similarity to the electrostatics. It 
has been shown that the director-mediated interaction is of a long range and possesses many 
other properties characteristic of the interaction between electric dipoles and quadrupoles 
UlElEllllElElElEllSlIinilll]. Particle trapping techniques [12] have been used to test this 
analogy and demonstrate experimentally the dipole-dipole [l3l[Tl], quadrupole-quadrupole 
flbl [16], and mixed disclination-dipole [17] pair interactions. Reorientation of the elastic 
dipoles was shown to be responsible for phase transitions between different 2-dimensional 
colloidal lattices on a nematic-air interface [18]. In the context of this analogy it is natural 
to expect that the Coulomb interaction, which is fundamental to the electrostatics, has an 
important implication in the physics of nematic emulsions, too. Recently we developed the 
electrostatic analogy in nematic emulsions to the level of charge and its density [IHl [20] . The 
director-mediated Coulomb-like interaction of two colloids was shown to be fully determined 
by vectors T^^^ and T^f^ of the transverse external torques (perpendicular to the unperturbed 
director at infinity) applied on the colloids [IHl [SO]- The scalar product — (r||^^ ■ r^,^'') 
plays the role of the product of two electrostatic charges in the 1/r interaction potential, 
and thus the two components of external torque play the role of two component elastic 
charge. Because of the difference between the scalar electrostatics and vector nematostatics, 
the elastic analogues of the surface charge density, charge, and higher multipole moments 
consist of two tensors (dyad). The multipole moments are naturally expressed via the 
elastic charge density which is determined by the two transverse director components on 
the surface imposing the director deformations. The interaction of the axially-symmetric 
sources, considered phenomenologically in [7] , obtains as particular case of the interaction of 
two correspondent multipole dyads. Small parameter of the theory is the ratio a/r = (colloid 
size/distance between colloids). For small a/r the theory provides all the tools available 
in the electrostatics, e.g., for solving different boundary problems that can occur in the 
nematostatics of anisotropic emulsions. In this paper we apply the nematostatics developed 
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in [ini [20] to the interaction between an elastic charge (dyad) and elastic dipole (dyad) with 
a wall (surface bounding the NLC) with different director alignments, which is the elastic 
counterpart of the well-known electrostatic problem solved by the method of images. In the 
next section we briefly introduce the colloidal nematostatics of Ref . fl9[ 120] and show that 
the integral form of the torque balance plays the role similar to that of Gauss' theorem in 
the electrostatics. Then the theory is applied to the colloid-surface interaction. A wall with 
homeotropic or planar director is shown to induce a repulsive force on the elastic dipole. 
The external torque, however, induces the elastic charge in this colloid and triggers switching 
to the repulsion. These results suggest that predictions of the colloidal nematosctatics 
can be tested by observing behavior of a single colloid at a sample surface which, in some 
situations, can be more robust than dealing with two colloids. 

II. ELASTIC CHARGE DENSITY REPRESENTATION OF THE COLLOIDAL 
NEMATOSTATICS 

A. Torque balance, Gauss' theorem, and elastic charge in 3 dimensions. 

The fundamental physical quantity of electric charge is purely phenomenological and 
must be postulated in the theory of elementary particles. In contrast, the nematostatics of 
the director field n allows for introduction of two different charges. Electrostatic potential 
is a scalar described by the linear Laplace (or Poisson) equation. It is the linearity that 
underlies the definition of the electric charge and its density as the source of electric field. 
At the same time, n is a vector field which reduces to a single variable, described by a 
linear equation (in the one constant approximation), only in 2 dimensions (2d). Owing to 
the linearity, the deformation source can be straightforwardly established: core of a point 
defect plays the role of a charge in 2d [211 ESI ESI El]- The independence of the integral, 
expressing the topological invariant, of the integration contour plays the role anlogous to 
Gauss' theorem in electrostatics, the invariant itself plays the role of a conserved charge, 
and the 2d nematostatics is similar to the 2d electrostatics with its logarithmic potential: 
disclinations of the same signs repel and those of the opposite signs attract each other. 

In 3d, however, the analogy between topological defects and charge is completely lost. In 
3d, the field n is described by highly nonlinear equations [21\ so that point defects, though 
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remain topological invariants, cannot be linearly connected with the distortions of n they 
induce [71 [TU]. Here the deformation sourse is the director distribution around the particle 
in its close vicinity of a size ~ a. We refer to such the deformation domain as particle though 
the distortion therein can be induced by surface of a real particle, by topological defects with 
zero total topological charge [TJ |T0], or by an external field dying out outside the domain 
area. Fig. 1. Consider 3-d director field n(r), uniform and parallel to the 2;-axis at infinity, 
Hoo = (0,0,1). At distances r ^ a, the small particle-induced perturbation rit of is 
transverse, t = x,y, and (in the one-constant approximations assumed in this paper) has 
the form 

Mr) = f + + 5 ^^^^^ ' + (1) 

It is natural to identify the coefficients with the subscript t in this expansion with the t-th 
component of elastic charge, elastic dipole, and elastic quadrupole, respectively. We seek 
the elastic analog of charge, following de Gennes' idea outlined in [2TJ. A transverse external 
torque = (r^,, Ty, 0) applied on a particle in the equilibrium is balanced by another, the 
elastic torque distributed over a surface S enclosing the particle, i.e., 

Tt = K Satpirpdprijdan^ + npdi3na)dSf3, (2) 
Js 

whrere K is the elastic constant, Satp is the absolute antisymmetric tensor, all indices but 
t run over 1,2,3, and summation over the repeated indices is implied. The integral in the 
r.h.s. does not depend on the choice of enclosing surface S, and the equality ^ reminds one 
Gauss' theorem with Ft in place of the electric charge. To further justify this connection 
one notices that integral ^ over a remote surface S vanishes for any term in the expansion 
Q but the first one. Substituting rit = qt/r in ^ and integrating over a large sphere gives 
Ty = 4:7iKq^, = -4:71 Kqy, or 

Thus, the tentative conclusion is that, in 3d, the role of Gauss' theorem and charge is 
played, respectively, by the balance of external and elastic torques (conservation of torque) 
and by the transverse components of the external torque exerted on the particle. This 
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is fully justified by calculating the Coulomb-like interaction in the elastic charge density 
representation developed in [121 120] • 




r=0 (elastic dipole) ^ r^O (elastic charge) 



FIG. 1: General deformation source ("partcile') in 3d. Inside the gray sphere the deformations can 
be very large (e.g., induced by point defects or strong anchoring of a real colloid). But outside the 
larger sphere the deformations are weak and linear which allows for the electrostatic anlogy The 
particle itself is of the dipolar type, but an external torque upon it can charge it, and it becomes 
an elastic charge. 



B. Dyads of elastic mutipoles and their interaction via the director field: the 
outline 

The results of Refs. [HI |20] instructive to our task here can be summarized as follows. 
Consider a deformation source with the director distributions rit given on the surface of 
enclosing sphere S with radius a. The quantity 

at(s) = nt(s)/a^ (4) 

plays the role of two component surface elastic charge density on the sphere 5*. Using natural 
analogy with the electrostatics, we define the two component elastic multipoles via the 
surface charge density as the following integrals over sphere S enclosing the particle: 
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9t = ^ / (^td^s, (5) 



s 



dt,a = ^ I ^tT^ad^s, (6) 



a3 



s 

2, 



Qi,a/3 = ^ (Tt{3Uaiyf3 - 6a(})d S, (7) 

Stt Js 

where u is the vector of a unit outer normal to S. With these definitions, we obtained the 
interaction of two particles with similar multipoles. Interaction of two "charged" particles 
with nonzero qt is Coulomb-like: 

^(1) . ^(2) . p{2)^ 

Ucoulom, = -^^K^-^ = - ^^^^ , (8) 

where we used relation (3), T^j^ = T^^^ — (r*^*^-noo) is the transverse component of the torques 
exerted upon the i-th particles, and R is the modulus of the separation vector R. The above 
connection between the elastic charge and external torque is thus fully justified. Eq. ([s]) 
shows that, depending on the sign of {T^f^ ■ T^f'), the elastic Coulomb interaction can be 
attractive or repulsive. In contrast to the electrostatics and 2d nematostatics, the charges 
with the same sign attract and with different signs repel each other ("parallel torques" 
attract whereas two " antiparallel torques" repel each other). Although the colloids must 
be anchored to the director, the Coulomb-like interaction does not directly depend on their 
specific shape and anchoring. Instead, the elastic charge is determined by the coefficients 
describing the torque exerted upon the colloid by a given type of external field. For instance, 
this can be the vector of permanent electric and magnetic dipole or electric and magnetic 
polarizability tensors of a given colloid. 

If the external torques are absent, the interaction energy is expressed solely in terms of 
particles' multipoles. The interaction between two "dipolar" particles is of the form 

where u = H/R is a unit vector along the separation direction. Eqs.([T]),(|4])-(|9]) (along with 
the quadrupole-quadrupole potential derived in [ElEO]) suggest the following interpretation. 
qt is the t-th component of the elastic charge and o"((s) is its surface density at point s on 
the sphere. The vector and tensor are the t-th dipole and quadrupole moments 
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determined in the standard way by the surface charge density at on the sphere. As (Tx and 
ay are separate sources, they determine not only the x and y director components outside 
the particle, Eq.Q, but also two independent tensors (dyad) for each multipole moment, 
i.e., qx and g^, d^^ and dy, Q,x and Qj^, and so on. 

III. COLLOID- WALL INTERACTION IN A NEMATIC LIQUID CRYSTAL 

The above formulas can be used to solve boundary problems similar to those of electro- 
statics. The simplest boundary problem is the interaction of an elastic multipole with a 
surface bounding the nematic sample and imposing planar or homeotropic director align- 
ment. Here we consider this problem for an elastic charge and dipole. 

A. Repulsion of an elastic charge from the wall 

Let us consider a single particle with a charge at a distance h from a plane surface of a 
NLC sample. We assume that the anchoring is strong, the director alignment in the sample 
far from the particle is homogeneous and parallel to the z-axis, rioo = (0,0, l),but the 
angle it makes to the surfaces is arbitrary. The charge qt can be induced by an external field 
exerting the torque T with the components Ty = AirKq^ and = —AirKqy. To justify the 
linearized theory, h is assumed to be large compared to the particle' size. As the director 
on the sample surface is fixed, the boundary condition is = 0, t = x, y. 

The problem can be solved using the mirror-image method. Let us place the image- 
particle with the charge q[ on the other side of a surface at distance h from it. Figs. 2, 3. For 
large h, distortions induced by the charge and its image are given by the sum (see Eq.([T])): 



where ri and r2 are the distances from a given point of the wall to the location of the charge 
and its image. As ri = r2, the boundary condition is satisfied for qt = — g^, t = x,y. The fact 



two opposite elastic charges repel each other, the elastic charge-wall interaction is repulsive. 
The repulsion force obtains from the interaction energy ([s]) of the torques T± and —T± by 
differentiating with respect to i? at i? = 2h, i.e.. 




(10) 



that the particle and its image are oppositely charged means that = -r^, Figs.2,3. As 
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The result depends on the direction of and thus on the surface tilt via the relation 



B. Repulsion of an elastic dipole from the wall 

Now consider the interaction between a wall and an elastic dipole represented by dyad 



nents. Symmetry makes some of them vanish. For ellipsoids with one of their axes along 
rioo = (0, 0, 1), the dipole dyad is diagonal: d^,. = (rf^., 0, 0), = (0, dy, 0) where ^ dy. In 
the case of an axially symmetric particle with symmetry planes passing through the sym- 
metry axis assumed to be along noo, d^ = dy = d (note that an axially symmetric particle 
without symmetry planes, such as a helicoid, is a chiral source which will be considered 
elsewhere). We restrict our consideration to this simple and practically important case of 
colloids. For instance, such are the so-called "topological dipoles", i.e., spherical particles 
with homeotropic boundary conditions with a companion hyperbolic hedgehog or disclina- 
tion ring [3 [TU] . 

From the general equation (|9]), the interaction energy of two axially symmetric dipolar 
particles with nonzero components d^^^ and d^"^^ obtains in the form 



where 6 is the angle the separation vector R, which, in our geometry, is along the surface 
normal, makes with the far homogeneous director rioo- This formula, up to the coefficient 3, 
reproduces the one obtained for the axially symmetric "topological dipoles" in [71 [TU]. 

Consider an elastic dipole dyad d^, t = x,y, aX a. distance h from the sample surface 
(wall). The z-direction is along the unperturbed homogeneous director noo = (0,0, 1). The 
general boundary condition on a surface with strong anchoring of any type is again n^ = 0, 
t = x,y. The image-dipole d^ is located at the distance h on the opposite side of the wall. 
The director field at point r of the surface is 



r^ = r-(r-noo). 



{dx,dy), Eq.(|6]). In general, each of the two vectors d^^ and d^, has three nonzero compo- 




(12) 
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„.M=3^*3^+3^i#, (13) 



where ri = and r2 are separation vectors between point r of the surface and the particle 
and its image. The boundary condition = gives two equations, i.e., 



(di • n) + (d; • r2) = 0, t = x,y. (14) 

We will consider the wall with the planar and homeotropic director alignment individually. 

a. Planar wall. The wall with a planar director alignment, Fig.4, coincides with the 
yz-plane a; = 0, while the x-axis is normal to the wall. Obviously, if ri = {—x,y,z), then 



r2 = {x,y,z). The two equations (14) then are solved by = d^ and d^, = — d^. Fig 
4. The interaction energy of the dipole dyad {dx,dy) and its image (d^.,— dj,) is readily 
calculated from eq.([9]) by substituting d^ = {d, 0, 0) and d^, = (0, d, 0). The force is obtained 
by differentiating this expression with respect to R at i? = 2h. This gives a repulsive force 
with the magnitude 

27d^ , , 

Fd,planar = AttK^^. (15) 

b. Homeotropic wall In the homeotropic geometry. Fig. 5, the uniform director and the 
2;-axis with the onset at the wall are normal to the wall which coincides with the xy-plane. 



Obviously, if ri = {x,y, —z) then V2 = {x,y,z). The two equations (14) then are solved by 
da; = — d^ and dy = — d^. Fig 5. The force is repulsive and has the magnitude 

9d'^ 

Fd,hom = ^^K—. (16) 



This force is 1.5 times weaker than F, 



d,planar- 



IV. CONCLUSION 



The nematostatics in 2 and 3 d is very different. The former is very similar to the 
2d electrostatics where disclination cores are in place of electric charges. The latter is 
similar to the electrostatics only in that its Green functions are Coulomb-like. In 3d the 
counterpart of the electric charge density is a dyad, the elastic charge can be induced only 
by an external torque whose components play the role of an elastic charge dyad. In this 
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3d colloidal nematostatics, the Coulomb-like interaction has the reverse sign. We described 
some implications of the colloidal nematostatics in 3d and showed that, in contrast to the 
electrostatics, the charges and dipoles are repelled from the wall. One interesting effect is 
that, applying the field-induced torque on a colloid, one can charge it, Fig.l. If the colloid 
is an elastic dipole, then by applying the external field one can switch the repulsion from 
the nematic surface from to 1/h? regime. Our results prompt the experimental tests of 
the interaction in nematic emulsions that, rather than dealing with a pair of particles, can 
deal with a single colloid at a wall. 
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FIG. 2: Elastic charge at a wall with fixed planar director alignment. Elastic charge qt induced 
by an external torque V and its image —qt induced by the image-torque — F . The director at the 
wall remains unperturbed and equal to noo. 
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FIG. 3: Elastic charge at a wall with fixed homeotropic director alignment. Elastic charge qt 
induced by an external torque F and its image —qt induced by the image-torque — F . The director 
at the wall remains unperturbed and equal to noo. 
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FIG. 4: Dyad of elastic dipole at a wall with fixed planar director alignment. The dipole (left) 
with da; = (d, 0, 0) and dy = (0, d, 0) and the image (right) with = (d, 0, 0) and d'y = (0, —d, 0) 
shown in two mutually perpendicular planes: a) xz-plane normal to the wall and b) planes x = —h 
(left) and x = h (right) parallel to the wall. 



13 




FIG. 5: Dyad of elastic dipole at a wall with fixed homeotropic director alignment. The dipole (left) 
with d^. = (d, 0, 0) and = (0, d, 0) and the image (right) with d^ = {—d, 0, 0) and d'y = (0, —d, 0) 
shown in the xz-plane normal to the wall. 
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